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Abstract 

The adsorption of oxygen on the Ag(100) is investigated by means of density func- 
tional techniques. Starting from a characterization of the clean silver surfaces oxygen 
adsorption in several modifications (molecularly, on-surface, sub-surface, Ag20) for 
varying coverage was studied. Besides structural parameters and adsorption ener- 
gies also work-function changes, vibrational frequencies and core level energies were 
calculated for a better characterization of the adsorption structures and an easier 
comparison to the rich experimental data. 
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1 Introduction 



Despite the progress realized in surface science during the last decades, there 
are still unsolved questions for processes as elementary as atomic adsorption. 
One particular example is the adsorption of oxygen on the Ag(OOl) surface. 

Already an early work-function study [1] indicates a non-trivial adsorption be- 
havior: while the work-function increases monotonically with oxygen coverage 
for the (110) surface of silver, the work-function change for the (001) surface 
depends not only on coverage, but also on temperature. For surface tempera- 
tures higher than about 320 K oxygen exposure results always in an increased 
work-function, whereas at lower temperatures only small oxygen doses lead 
to a positive work-function change. For higher exposure the work-function de- 
creases upon oxygen adsorption. Later this "phase transition" between a high- 
and a low-temperature structure could be reproduced and the different phases 
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characterized by high resolution electron energy loss spectroscopy (HREELS) 
and low-energy electron diffraction (LEED) [2]. Very recently these results 
were confirmed by a combined LEED, HREELS, x-ray photoemission spec- 
troscopy (XPS) and x-ray photo-electron diffraction (XPD) study [3]. Upon 
oxygen exposure several distinct species can be prepared (comp. table f): 

(f ) Adsorption at low temperature (T < 150K) leads to a molecular species, 
which is stable up to 200 K and has been well characterized by XPD [3], 
HREELS [7] and thermal desorption spectroscopy (TDS)[8]. According to 
a scanning tunnelling microscopy (STM) study, the molecules adsorb in 
hollow sites, forming a c(2 x 4) structure[6]. 

(2) After dissociation atomic oxygen adsorbs at the surface with an O Is 
binding energy of 530.3 eV and a vibrational frequency of about 30 meV. 
At around 190 K a sharp c(2 x 2) pattern appears and the frequency shifts 
up to 36 meV. An additional loss at 131 meV has been assigned recently 
on the basis of isotopically labelled oxygen to an electronic transition [5]. 
Based on XPD a missing-row p(2a/2 x y/2) reconstruction of the substrate 
was proposed (with the oxygen atoms forming a c(2 x 2) arrangement) [3]. 

(3) Heating of this structure above ~ 320 K lifts the reconstruction, the O 
Is binding energy changes from 530.3 eV to 528.3 eV and the vibrational 
frequency decreases again from 34-37 meV to 28-31 meV [3]. 

(4) Around the transition temperatures also an additional oxygen moiety 
appears with a binding energy of 530.9 eV [3], which vanishes already 
at 350 K and could have been shadowed for lower temperatures by the 
dominant peak at 530.3 K. This species has been assigned to sub-surface 
oxygen. 

Interestingly only the higher temperature species characterized by an O Is 
binding energy at 528.3 eV is active towards CO and C 2 H 4 oxidation[3]. Based 
on the structural models from XPD, a combined ab-initio and ultraviolet pho- 
toemission spectroscopy study was devoted to the investigation of electronic 
surface states [9]. In a recent STM experiment three different species of ad- 
sorbed oxygen (T a d s = 150 K) are reported for low coverage (0o ~0.1 ML) and 
assigned to hollow, bridge and on-top adsorbed atoms [10]. Finally, last year 
two density functional studies appeared investigating the oxygen adsorption 
on the Ag(001) [11,12]. Both studies are in a good agreement with our study. 
Also the adsorption on Ag(lll) surface was investigated recently via DFT 
techniques [13,14]. 

In the present study we use density functional theory (DFT) calculations 
to gain further insight into this complicated adsorption behavior. First we 
present our results for the clean silver surface. Then we discuss adsorption of 
oxygen in high symmetry (on-surface) sites of the (001) surface at different O- 
coverage. Subsequently results for molecular adsorption, the proposed missing- 
row reconstruction and the silver oxide surface are presented. The paper closes 
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Table 1 

Experimental characterization of the different oxygen adsorption phases on Ag(001): 
O Is binding energies E b {0 Is), Ag-0 stretching frequencies his, LEED pattern, 
work-function change A<3> and estimated height of O above the surface ho- The 
results for atomic oxygen are taken from Refs. [1,2,3,4,5], those for molecular oxygen 
(T< 200 K) from Refs. [6,7]. 

with a discussion of our results in the light of the experimental data. 



2 Methodology 

The calculations are performed using the Vienna ab-initio simulation pack- 
age VASP [15,16] which is a DFT code, operating in a plane-wave basis set. 
The electron-ion interaction is described using the projector-augmented-wave 
(PAW) method [17,18] with plane waves up to an energy of E cut = 250 eV. For 
exchange and correlation the functional proposed by Perdew and Zunger [19] 
is used, adding (non-local) generalized gradient corrections (GGA) following 
Perdew et al. [20]. If not mentioned differently we have used a (12 x 12 x 1) 
k-points mesh for the integration over the Brillouin zone for the p(lxl) cell, 
(6x6x1) for c(2x2), (4x4x1) for p(2x2) and a (3x3x1) mesh for a p(3x3) 
cell. 

The system itself is modelled by six layers of Ag separated by 14 A of vacuum 
with oxygen adsorbed on one side of the slab. The upper two layers of the Ag 
surface are allowed to relax, the remaining atoms are fixed at their ideal bulk 
positions (acGA=4.16 A, a ex p=4.08 A). 

The free molecule is characterized by a calculated stretch frequency of 1561 
cm" 1 at an equilibrium bond length of 1.23 A. The corresponding experimental 
values are 1580 cm _1 and 1.21 A[21]. The problematic description of the O2 
binding energy (Eq^ — 6.25 eV, E^f =5.23 eV) stems mainly from the error in 
the energy of the free atoms, where high density gradients make an accurate 
description more difficult. So the reported energies (with respect to the calcu- 
lated binding energy of molecular oxygen) should be of much higher accuracy. 



3 



This is for example reflected by the calculated formation energy for Ag20 of 
0.39 eV compared to the literature value of 0.32 eV[22]. 

Adsorbate-substrate stretching frequencies have been obtained from finite dis- 
placements of the oxygen atoms along all three Cartesian directions (A = 
±0.04 A). In order to increase the accuracy we performed the calculations with 
a harder pseudo-potential requiring an energy cut-off of E cut = 400 eV. Diffu- 
sion barriers were determined with the nudged elastic band method [23,24]. 

Core-level shifts were calculated in a (2y/2 x \/2) cell sampled by a (3 x 6 x 1) 
k-points mesh. Further details concerning the calculation of the core levels 
(final state effects) are given in Ref. [25]. 



3 Results 

3.1 Clean Ag (001) 

The clean Ag(001) surface is well characterized, experimentally as well as 
theoretically [11,12,26,27,28,29,30,31]. Calculations have been done for slabs 
with 4 to 8 layers, to check for errors due to the finite thickness of the Ag 
film. Both surfaces of the slab were relaxed. Table 2 compiles our results for 
inter-layer relaxation, surface energy and work-function together with values 
from the literature. We find that structural, as well as electronic parameters 
agree well with earlier results. A slight inward relaxation of the surface plane 
is compensated by outward relaxations of the lower layers of comparable size. 
The surface energy is underestimated by about 30 % compared to earlier LDA 
calculations and experiments, a property of the GGA that has already been 
discussed earlier [32]. However, since it is well established that the GGA leads 
to a more reliable description of the adsorption process we did our calculations 
with the inclusion of gradient corrections. 

Concerning the slab thickness dependence, results for the six layer slab, which 
were used for the adsorption studies are converg ed within 0.1% (0.02 A) for 
the inter- layer relaxation, within 30 mJ/m 2 for surface energies and the work- 
function varies only by 60 meV when going from six to eight layers. 

3.2 Adsorption in high symmetry on-surface sites 

As a first step towards a deeper understanding of the adsorption process, we 
put oxygen at varying coverage (6 = |, |, |, |, 1 ML) in the high symme- 
try positions (top, bridge, hollow) of the slab, optimized the geometry and 
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Table 2 

Inter-layer relaxation (Ac?i2, Acfe, ACZ34) in %, surface energy (a) in J/m 2 and 
work-function (<£) in eV for a clean Ag(001) surface as a function of slab thickness. 
Earlier theoretical and experimental results are taken from Refs. [27] (a), [28] (b), 
[11] (c), [12] (d) and [29] (e), [30] (f), [31] (g). 

calculated frequencies and work-function changes. For the characterization of 
the energy balance during adsorption we chose two slightly distinct quantities. 
First, the adsorption energy per adsorbed atom (E^s) with respect to half an 
oxygen molecule defined as 

N 

Ea.ds = (E tot - -E c l can — • E 02 )/N , (1) 



with the total energy of the adsorption system (containing iVo oxygen atoms) 
E tot , the energy of the clean surface (-Edean) and the energy of a free 2 
molecule Eq 2 . Following this definition, E ads decreases (i.e. becomes more 
exothermic) with increasing Ag-0 bond strength. Secondly, we use a (gen- 
eralized) surface energy per area A ( especially for the comparison with the 
reconstructed missing-row structure, discussed below) 

No 

c = (-^tot - N Ag ■ E Ag — • Eq 2 )/A - <7 c i ean , (2) 

where A^ g and E Ag are the number of silver atoms per cell and the total energy 
of a bulk silver atom, respectively. The subtraction of the surface energy of a 
clean silver surface <7 c i ea n accounts for the second (uncovered) surface of the 
slab. This surface energy describes the energy cost for creating a surface minus 
the energy gain due to oxygen adsorption. 

The results for the various adsorption sites at different coverage are compiled 
in table 3; fig.l gives a graphical summary of the adsorption energy versus 
adsorption height. At every coverage, adsorption in hollow sites is favored while 
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on-top and bridge adsorption is nearly always endothermic. Only at a coverage 
of = 0.25 ML adsorption in bridge sites becomes thermoneutral, but remains 
still ~ 0.8 eV less favorable compared to adsorption in the hollow sites. All 
these findings for the adsorption energetics are in nice agreement with two 
recent DFT studies [11,12]. This strong preference of hollow adsorption makes 
the interpretation of the STM images in Ref. [10] in terms of bridge- and 
top-adsorbed meta-stable oxygen species very improbable, unless there were 
huge barriers between the adsorption sites. This point will be addressed in the 
following section. 

At full coverage we performed additionally spin-polarized calculations. For 
atop adsorption a substantial magnetic moment of the O atom of 1.4 /xb per- 
sists, leading to a decrease in the adsorption energy by -0.4 eV at an unchanged 
geometry. In the higher coordinated sites the effect of spin polarization is neg- 
ligible, due to the stronger interaction with the substrate. At lower coverage 
this effect is expected therefore to be even smaller due to stronger Ag-0 in- 
teraction, justifying the neglect of spin polarization in further calculations. 

The inward relaxation of the first substrate layer, which was observed for the 
clean surface, persists up to a coverage of 6 = | ML; for higher coverage 
it is gradually turned into an expansion. For the second inter-layer spacing 
the sign of the relaxation changes already at G = \ ML. The adsorption 
heights for the three high symmetry adsorption sites vary such, that the Ag-0 
bond length remains always between 1.9 and 2.3 A. This results in adsorption 
heights of around 2.0 A, 1.4 A and 0.7 A for top, bridge and hollow adsorp- 
tion, respectively. The correlation between adsorption energy and adsorption 
height reveals a very peculiar character of O-adsorption in the hollow position 
(see fig. 1). Whereas for bridge- and top-adsorbed O the adsorption-height is 
almost independent of the coverage and the adsorption energy becomes more 
endothermic as the coverage increases, hollow-adsorbed O sinks deeper into the 
surface with increasing coverage - at full coverage the stable adsorption sites 
are located even below the surface-plane (/ie=i — —0.33 A ). This reduction 
of the adsorption-height is accompanied by a strong outward relaxation of the 
top Ag-layer (Arf 12 =31%). This expanded adsorption structure may be inter- 
preted as the formation of an silveroxide layer which is weakly bonded to the 
Ag-substrate. A similar behavior is found already for = | ML (p(2x2)30). 
In this structure two inequivalent adsorption sites exist for oxygen. Two out 
of three oxygen atoms are equivalent due to symmetry and stay above the 
surface at a height of h = 0.26 A, while the third oxygen atom sinks below 
the first layer (h = -0.09 A ). 

The work-function change A$ for top and bridge adsorbed oxygen increases 
monotonically with coverage. For the hollow site this trend holds only up to 
| ML, from where on A$ decreases again, due to the aforementioned negative 
adsorption heights at 0.75 ML and 1 ML. 
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The vibrational frequency perpendicular to the surface exhibits for all three 
adsorption sites an increasing trend with coverage up to around 0.5 ML. From 
there on the vibrational modes become softer again due to the pronounced 
weakening of the Ag-0 bond at high coverage. 
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Fig. 1. Adsorption energy {E a( is) versus height of the adsorbed O atom for various 
adsorption sites and coverage. The second data set for full coverage (denoted by 
an asterisk) has been obtained from a spin polarized calculation. For subsurface 
adsorption also energies for incorporation in bulk Ag are included (empty circles). 



3.3 Transition between high- symmetry sites - surface diffusion 



In a recent low temperature STM study [10] three different oxygen species 
have been imaged for low oxygen coverage. These have been assigned to hol- 
low adsorbed oxygen atoms and two meta-stable species residing in bridge and 
on-top sites, respectively. Considering the large differences in the adsorption 
energies for these three sites, this could only be possible if there were significant 
barriers separating the meta-stable sites from the hollow position. In order to 
clarify this point we performed nudged-elastic-band (NEB) [23,24] calculations 
to determine the diffusion barriers between neighboring high symmetry sites. 
Fig. 2 shows the variation of the potential energy for the movement of an oxy- 
gen atom from a bridge-site via hollow towards the on-top position, following 
the pathway shown in the inset. For this calculation (performed in a p(2x2) 
cell) again the upper two silver layers were allowed to relax. 
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Table 3. Results for oxygen adsorbed in high symmetry sites (t - top, b - bridge, h - hollow) on Ag(001) for coverages varying between 
= 0.11 ML and = 1 ML. The reported quantities are: a - surface energy, E a d s - adsorption energy, h - height of O with respect to 
the average surface plane atoms, Ad\2, Ac?23 - change of the average inter-layer spacing, 61,62 - buckling of I s * and 2 nd layer, dAg-o - 
minimal Ag-0 distance, A$ - oxygen induced work-function change, v - oxygen frequency perpendicular to the surface. 
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Fig. 2. Adsorption energy E a( i s along a pathway connecting the adsorption sites 
bridge - hollow - top calculated in a 2 x 2 cell (octahedral - hollow - tetrahedral 
sites for the pathway leading into subsurface positions; compare sketch). Calculated 
points are indicated by crosses, the connecting lines are spline fits, based on forces. 

The results show that bridge and on-top sites are only saddle-points of the 
potential energy surface for this coverage. Hence on the basis of these findings 
the interpretation given in the STM work becomes very improbable. However, 
the calculation reveals another possible (3-fold) adsorption site between hollow 
and on-top, reflected by a dip along the pathway (arrow in the graph on fig. 2). 
For a further investigation we put the oxygen atom artificially at that position 
and performed again a structural optimization. Our calculations did not lead 
to any local minimum. 



3.4 Sub-surface adsorption 



We have also considered the sub-surface adsorption of oxygen, since such 
species are believed to be important for the selectivity of silver in catalyz- 
ing ethylene epoxidation. [33] We investigated the sub-surface adsorption for 
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coverages between a quarter and full monolayer (p(2x2), c(2x2) and p(lxl) 
structures) in the octahedral and tetrahedral sites, the results of which are 
compiled in table 4 and fig. 1. For comparison, we included in fig. 1 also 
the energies for oxygen incorporated in the tetrahedral and octahedral sites 
of bulk Ag (in a 3 x 3 x 3 cell). Sub-surface adsorption is an endothermic 
process compared to molecular O2 at all investigated coverages (see fig. 1). 
However, there are interesting trends of the adsorption energy in coverage 
and in comparison to the dilute solution of oxygen in the bulk. In contrast to 
the on-surface adsorption the variation of the adsorption energy with cover- 
age is non-monotonous, the smallest (least endothermic) energy is calculated 
for = 0.5 ML, i.e. at a composition of the surface layer corresponding to 
stoichiometric Ag20. And the dilute solution of oxygen in bulk silver is an 
extremum for both adsorption sites: for the octahedral site it is more stable 
than the high coverage subsurface case, whereas the tetrahedral interstitial is 
even less stable than the quarter coverage subsurface position. 

This peculiar behavior can be understood when considering the geometries of 
the two interstitial sites. Placing an oxygen atom into an ideal (unrelaxed) 
Ag bulk leads to Ag-0 bond lengths of 1.80 A for the tetrahedral and 2.08 
A for the octahedral site. So compared to the typical Ag-0 bond lengths of 
2.0-2.2 A (compare table 3) the octahedral site can easily accommodate an 
oxygen atom, whereas the tetrahedral site is too small. This explains already 
the strong preference for the octahedral site in the dilute solution of oxygen in 
the bulk. At the surface the incorporation of oxygen induces pronounced ge- 
ometrical changes, partially just for sterical reasons, so that the site becomes 
large enough to accommodate the additional atom; partially because the for- 
mation of strong Ag-0 bonds weakens the Ag-Ag bonds and leads hence to 
an outwards relaxation of the first layer. 

For the larger octahedral site, the second contribution dominates: already at 
the lowest coverage (p(2x2)), the outward relaxation together with the up- 
shift of the immediate silver neighbor in the first layer (reflected by a large 
buckling of 0.42 A), expands the octahedral site such that two of the six Ag-0 
bonds increase beyond 2.3 A and are significantly weakened. This results in 
the weakest subsurface adsorption energy of this study. With increasing cover- 
age, the outward relaxation increases dramatically, so that a slightly different 
adsorption geometry between the first and second layer becomes available. 
In this position the Ag-bond to the third layer is now completely absent 
(d l Ag _ > 3 A), but for that the remaining 5 bonds are stronger and make this 
position lower in energy. 

The smaller tetrahedral site on the other hand has to rely on the (local) 
expansion of the substrate. At low coverage the two silver neighbors in the 
first layer have to shift up significantly (leading to a buckling of 0.64 A), so 
that the interstitial becomes large enough. Such a large distortion costs energy 



10 





pU 


xlj 


c(2x 


2) 


p(2x 


2) 




o — 


1 1\/TT 


e = 0.5 ML 


9 = 0.25 ML 




oct 


tet 


oct 


tet 


oct 


tet 


a [J i m \ 


1.00 


i /in 
1.4y 


1.24 


1.12 


1.13 


1.08 


Kds [eVJ 


0.46 


u.o r 


n a a 

U.44 


U.ol 


n 70. 
u. < y 


n fio 
U.Oy 




£fl A 

0U.4 


04. z 


oO.O 


A A Cl 

44. y 


Q 1 
O. ( 


10. U 




n 


n 


O.oo 


n 
U 


n A o 
U.4Z 


n /i 
0.64 


A<7oo [%1 


n q 

u.y 


— n q 

u.y 


1.4 


0.6 


3.0 


-0.5 


/i„ r A i 


u 


n 
u 











0.12 


d Ag-0 i A \ 


z.lo 


Z.I / 


2.21 


2.15 


2.38 


2.10 


a™ r A i 


Z.oU 




2.27 




2.23 




d^-o [A] 


3.07 


2.19 


3.02 


2.20 


2.32 


2.18 


A$ [eV] 


0.15 


0.37 


0.15 


0.37 


0.27 


0.15 


[meV] 


34.2 


37.4 


36.0 


48.8 


29.5 


49.6 



Table 4 

Results for oxygen adsorbed in the octahedral (oct) and tetrahedral (tet) sub-surface 
sites at coverages of 1 ML (p(lxl)), 0.5 ML (c(2x2)) and 0.25 ML (p(2x2)). The 
reported quantities are: a - surface energy,-E ac ; s - adsorption energy, AeZr2, Ac?23,- 
change of the average inter-layer spacing, 61,62 - buckling of I s * and 2 nd layer, d Ag _ Q , 
d r j ig _o-, d l A g -o ' shortest Ag-0 distance to the upper, middle and lower Ag layer, A<1> 
- the work function change compared to the clean Ag(001) surface and the oxygen 
frequency perpendicular to the surface v. 

and makes the tetrahedral site unfavorable at low coverage. With increasing 
coverage, the weakening of the Ag-Ag bonds opens up the first interlayer 
spacing, so that adsorption in the tetrahedral site becomes not only possible, 
but even more favorable than adsorption in the octahedral site. 

The vibrational frequencies perpendicular to the surface for the tetrahedral 
and octahedral site are remarkably similar to those for the corresponding on- 
surface sites, bridge and hollow, respectively. 

Similar to the on-surface case we calculated the potential energy along a path- 
way connecting the subsurface sites with the on-surface hollow position in a 
p(2 x 2) cell. The monotonic decrease in energy towards the hollow position 
demonstrates that isolated sub-surface sites are highly unstable at this cover- 
age and can only be stabilized by co-adsorbed on-surface oxygen, preventing 
the migration back to the surface. 
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3.5 Missing-row reconstruction 



For the low-temperature phase (below 320 K, compare table 1), a missing- 
row reconstructed phase with a p(2v^2 x y/2) periodicity has been proposed 
by Rocca et al. on the basis of XPD experiments [3] . The still rather regular 
position of the oxygen atoms should give rise to the observed c(2x2) LEED 
pattern [2]. 
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Fig. 3. Results for the (2v2 x v2) missing-row reconstruction with two oxygen 
atoms in the pseudothreefold hollows. The sketch shows a top and side view of the 
optimized structure (all distances in A). The reported quantities in the table are: 
a - surface energy, E ac [ s - adsorption energy with respect to a clean reconstructed 
surface, h - height of O with respect to the average surface plane, Adi2, Ag?23 
- change of the average inter-layer spacing, 61,62 - buckling of I s * and 2 nd layer, 
dAg-o - minimal Ag-0 distance, A<I> - oxygen induced work-function change, v - 
vibrational frequency perpendicular to the surface. 

Starting from the proposed geometry we optimized the structure for the missing- 
row reconstruction (see fig. 3). However, in the structure proposed on the basis 
of the XPD experiments the oxygen atoms are 0.15 A below the average of 
the surface plane and laterally shifted by 0.36 A towards the missing rows [3], 
at variance to our results. 



In the relaxed configuration the oxygen atoms are located at a height of about 
0.13 A above the average surface plane in the pseudo-hollow sites along the 
missing-rows; the under-coordinated Ag atoms move by ~ 0.2 A towards the 
missing-rows, leading to a structure comparable to what is observed for O on 
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Table 5 

Several combinations of sub-surface oxygen(s) in the missing-row (2\/2 x \/2) cell 
in addition to two oxygens in the pseudo- hollows (comp. fig. 3). The first column 
describes the position of the sub-surface atoms by the atom number (according to 
fig. 3) of the nearest on-surface atoms. Tetrahedral sites (tet) are between first and 
second layer, octahedral sites (oct) in the second layer. Compiled is the resulting 
total coverage 0, the adsorption energy with respect to a clean reconstructed surface 
E a d s and the generalized surface energy a. 

Cu(OOl) [34,35]. Very similar results were obtained in the recent theoretical 
studies by Cipriani et al. [11] and Wang et al. [12]. 

Besides the differences in structure, the reconstructed phase also cannot ac- 
count for other properties observed for the low temperature phase. Experimen- 
tally this phase should be characterized by an increased vibrational frequency 
around 34-37 meV [3,5] and a negative work function change [1], compared to 
the calculated values of 31 meV for the frequency and a work function change 
of A$=+0.5 eV. 

From an energetic point of view (surface energy), this reconstructed structure 
is very similar to the simple c(2x2) structure. In order to minimize any errors 
due to differences in the setup (k-points, grids for Fourier transformations) 
we recalculated both structures (the c(2x2) and the missing row structure) 
in a (2\/2 x \/2) cell with a denser k-point mesh (5x10x1). Using this setup, 
the surface energy (calculated according to equation 2) for the reconstructed 
structure is even slightly lower in energy (Act =50 mJ/m 2 ), so that within the 
accuracy of our calculation both structures are essentially degenerate. 

Hence, although oxygen adsorption leads indeed to a strong weakening of the 
Ag-Ag bonds, so that structures such as the missing-row reconstruction be- 
come energetically possible, this particular structure is not consistent with the 



13 



experimental findings. A way out of this dilemma would be a more compli- 
cated structure with additional sub-surface oxygen, as it was already argued 
in Refs. [3,5]. Therefore we tested several combinations of occupied octahedral 
and tetrahedral sub-surface sites together with the experimentally proposed 
pseudo-hollow site along the missing- rows (see table 5). However, neither of 
these structures can either account satisfactorily for the observed experimen- 
tal structure nor is energetically more favored than the simple missing row 
reconstruction. 



3. 6 Molecular Adsorption 



Adsorption below 150 K leads to molecular adsorption of oxygen [36,7]. At 
low coverage two distinct molecular species could be identified via vibrational 
spectroscopy [36]. However, since already for coverages as low as 0.15 ML one 
species starts to dominate significantly, we will concentrate in the following 
only on this moiety, which saturates finally in a c(2 x 4) structure, as ob- 
served by STM [6]. By comparison with tight binding molecular dynamics 
simulations with these STM results the four-fold hollow site was proposed 
as adsorption site with the molecular axis oriented towards the neighboring 
bridge sites. The analysis of thermal desorption spectra [8] lead to an esti- 
mate for the adsorption energy of 0.4 eV. This molecular species is further 
characterized via two dipole active modes at 84 meV and 30 meV [7] and an 
O Is core binding energy of around 532.2 eV [3] (comp. table 1). We tested 
several molecular adsorption configurations. The most stable position of the 
molecule is over a hollow site, with the molecular axes oriented towards the 
bridge positions, in agreement with the proposed structure from Ref. [6], fig. 
4. In the c(2 x 4) structure (i.e. 0.25 molecules per substrate atom) every 
silver atom binds to only one molecule, so that a quite strong binding is pos- 
sible. The results of our calculation are compiled in table 6. The molecule is 
stretched to a bond-length of 1.43A, situated 1.47A above the surface plane. 
Similar to atomic adsorption, the first interlayer spacing is reduced by —0.7%, 
whereas the second is slightly expanded (+1.1%). The magnetic moment of 
the molecule vanishes completely and the work-function of the surface is in- 
creased by 1.83 eV. The calculated adsorption energy (-0.68 eV per molecule) 
is significantly higher than the experimental estimate of -0.4 eV[8], reflecting 
the well known overestimation of adsorption energies even with the used gra- 
dient corrected exchange correlation functionals [37]. Also the frequencies are 
overestimated by about 15%. A possible reason for this could be the neglect 
of substrate movement when calculating the frequencies, which could be for a 
"soft" substrate with a low melting temperature like silver quite crucial. 
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Table 6 

Results for molecular oxygen adsorbed in the hollow sites of Ag(001) (with the 
molecular axis pointing towards the bridge sites) in a c(2x4) structure. The reported 
quantities are: a - surface energy, E a ^ s - adsorption energy per molecule, molecular 
bond length (d) and height above the surface plane (h), A.d\2, Ac?23 " change of the 
average inter-layer spacing, 61,62 - buckling of 1 st and 2 nd layer and the change of 
the work- function A<&. The frequencies (u) correspond to intramolecular stretch, 
molecule-surface vibration, frustrated rotation perpendicular and parallel to the 
surface, and the two frustrated translations parallel to the surface. 

3.7 Silver-oxide - Ag 2 

Since it is nowadays well established that on some metals the reactive surface 
is the oxidized one [38,39], we performed also a limited set of calculations for 
the (001) surface of Ag 2 0, which is the stable oxide at room-temperature (see 
Fig. 4). The lattice constant as determined by the calculation is a GGA =4.83 
A(a exp = 4.72A[40]). In the DFT silveroxide exhibits a vanishing energy gap 
and becomes hence metal-like [41] in contradiction to experimental findings, 
according to which Ag 2 should be a semiconductor. However, we hope that 
these shortcomings in the electronic description do not have a major influence 
on the results presented in the following. 

For the modelling of the surface we chose a rather small super-cell with two 
silver atoms per layer and a total of six Ag-layers as sketched (after relaxation) 
in Fig. 4. With this setup we could model at the same time a silver terminated 
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Fig. 4. Structure of silver oxide (Ag20) in the bulk and of the (asymmetrically 
terminated) (001) oriented surface slab, as it was used in the calculations. 
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Table 7 

Results for the (001) surface of silver oxide (Ag 2 0): Change of the inter-layer 
spacing (Ad) with respect to the ideal (bulk) Ag-0 layer distance and the work- 
function <E>. 

surface on one end of the slab and an oxygen terminated one on the other side. 
We fixed the innermost two layers at their ideal bulk position and optimized 
all other layers. The results are compiled in table 7. Further results for this 
oxidic surface are given in the subsequent section. 

3.8 Core Level Binding Energies 

The phase transition between the low-temperature and high temperature phase 
has been characterized mainly by core-level spectroscopy (XPS) [3] (comp. ta- 
ble 1). In addition to the O Is binding energies also the Ag 3d$/2 core levels 
have been measured as function of coverage and temperature. Hence, we calcu- 
lated also for some structures the binding energies of the core-levels. In these 
calculations vanishing screening was assumed, so that all electrons can simul- 
taneously relax into the ground state, after the core electron was removed 
(final state effects) [25]. In order to minimize errors due to different unit-cells, 
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all calculations were performed within a the p(\/2 x 2y / 2) cell sampled by a 
grid of 6 x 3 k-points. Since the absolute values from the calculations are not 
well defined, only energy differences can be compared. For silver the bulk core 
level (calculated with in the same super cell) is an appropriate reference, for 
oxygen we chose hollow adsorption at low coverage (6 = 0.25 ML). 

3.8.1 Ag3d 

Experimentally, for the Ag 3d 5 / 2 core level only two peaks could be resolved 
independent of temperature and oxygen coverage: One peak at 368.0 eV for 
the clean surface component and a second at 367.6 eV for silver atoms in 
contact with oxygen. [3] While for the low temperature phase both peaks 
are almost of equal intensity, the surface component clearly dominates at high 
temperatures. Fig. 5 summarizes the results of our calculations. On the clean 
surface the Ag 3d 5 / 2 energy is shifted by approximately 240 meV with respect 
to the bulk Ag core-level. With increasing oxygen coverage the binding energy 
shifts down, almost linearly with the number of oxygen neighbors, as it was 
observed already for oxygen on Ru(0001) and Rh(lll) [42]. Since already 
a single oxygen neighbor induces a shift (with respect to the clean surface) 
of about 0.25 eV, the experimentally observed shift of —0.4 eV implies that 
there are no silver atoms with more than two oxygen neighbors. This is also 
in line with the experimental interpretation of a maximum coverage of 0.4 
ML [3]. Slightly different is the situation for the missing row reconstruction. 
Here the reduced number of silver nearest neighbors shifts the core-level up by 
approximately 150 meV, thus compensating partially for the oxygen induced 
shift (comp. fig. 5, full symbols). 

3.8.2 Ols 

The experimental results for the oxygen Is core-levels from Rocca et al. [3] are 
more complicated (comp. table 1, fig. 6): At lowest temperature the molecular 
species can be identified via its binding energy around 532.0 eV. Dissociation 
of the molecules leads to a prominent peak at 530.3 eV up to a temperature of 
around 300K, where this peak starts to vanish. Prior to disappearing another 
small peak appears at approximately 530.9 eV, which could have been there 
already before, but buried by the main peak at 530.3 eV. This peak vanishes 
also at 350 K. Above 300K the high-temperature species gives rise to a peak 
at around 528.3 eV. 

Table 8 summarizes the results for various structures and adsorption sites. All 
energies are given with respect to the O Is binding energy of hollow-adsorbed 
oxygen at low coverage (O = 0.25 ML). The core-level is mainly determined by 
the coordination of the adsorption site and changes hardly with coverage, as 
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Fig. 5. Ag 3d surface core-level shifts as a function of the number of oxygen neigh- 
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Fig. 6. Experimentally determined O Is binding energies for varying temperature, 
taken from Ref.[3] (left panel). Calculated core-level binding energies for different 
coverages and adsorption geometries (right panel); numbers indicate the coverage 
in mono-layer, dashed lines refer to oxygen in bulk positions. 

can be seen from the comparison of the core-level for hollow adsorbed oxygen 
at = 0.25 and 0.50 ML. Only for O = 1.00 ML the binding energy is larger, 
since the special adsorption geometry slightly below the surface layer for this 
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Table 8 

Oxygen Is core level energies for oxygen adsorbed for various coverage, sites (h- 
hollow, b-bridge,t-top,tet - tetrahedral subsurface site, oct - octahedral subsurface 
site) and structures. All energies are given with respect to the binding energy of the 
hollow adsorbed species at low coverage. 



coverage leads to a higher "effective coordination number". Also the values 
for the missing row reconstructed surface have slightly higher binding energies 
compared to hollow adsorption. Again the low adsorption height increases the 
effective coordination and hence also the O Is binding energy. The subsurface 
sites continue this trend: the oxygen binding energy of the lower coordinated 
tetrahedral site is 0.2 eV below that of the hollow adsorbed species, for the 
higher coordinated octahedral site the binding energy increases even more. 
The highest binding energy is calculated for the molecular state: in addition 
to the high-coordinated adsorption site there is a substantial oxygen-oxygen 
bond. The lower coordinated bridge and on-top adsorbed atoms continue the 
trend at the other end of the spectrum. 
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3.9 Optimal adsorption sites with increasing coverage: mixed surface and sub- 
surface adsorption 



To complete our study we investigated several more adsorption structures (in 
total more than 50 different structures) with mixed on-surface and sub-surface 
oxygen for total oxygen coverage between = 0.11 ML and = 2 ML. The 
adsorption energy for a selection of adsorption structures as a function of 
oxygen coverage is plotted in fig. 7. Up to half monolayer coverage the hollow 
site is the most favorable adsorption position; at = 0.5 ML only the missing- 
row reconstruction is degenerate with the c(2 x 2) structure (comp. Sec. 3.5). 
At higher coverage combinations of sub-surface and on-surface sites become 
more stable: at = 0.75 ML a combination of two hollow and one octahedral 
sites in a p(2 x 2) cell, at full coverage half of the oxygen is located below 
the surface in octahedral sites and half of it on the surface in the hollow sites. 
Interestingly for adsorption in isolated sub-surface sites, the tetrahedral site 
is more stable (see fig. 2). This is also true in combination with 1 ML of 
on-surface oxygen (see fig. 7), but the octahedral site becomes clearly more 
favorable for higher sub-surface coverage. 




_ "J Q I 

0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00 

coverage [ML] 

Fig. 7. Adsorption energy E a d s for the structures at the coverage from 0.25 ML to 
2.0 ML calculated in a p(2 x 2) cell. The second favorite structure is indicated by 
squares. 
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4 Summary and Discussion 



4-1 Surface energy 
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Fig. 8. Adsorption energy per area with respect to the chemical potential. Only the 
most favorite structures for each coverage (comp. fig. 7) are shown). Additionally 
included is the half coverage missing row structure with a dashed line. ^fi(02) = eV 
corresponds to free molecular oxygen and /i(Ag20) = -0.39 eV to silver-oxide. 

In order to facilitate a direct comparison of all structures with a different 
number of Ag atoms in the surface (missing row reconstruction, silver oxide 
Ag 2 0, . . .), we use the (surface) grand potential, i.e. the Legendre transform of 
the generalized surface energy (Eqn. 2) with respect to the number of oxygen 
atoms 



( \ N ° 



(3) 



where \io stands for the chemical potential of oxygen (with [loiP^) = eV). 
Via this transformation the dependence on the number of (oxygen) particles 
is transformed into a dependence on the chemical potential. So every struc- 
ture is represented by a line, where the intercept at the ordinate equals the 
surface energy and the slope increases with coverage. In Fig. 8 this quantity 
is plotted. It describes the stability of different Ag-0 structures as a function 
of varying oxygen chemical potential. For fixed temperature (partial pressure) 
the chemical potential scale is equivalent to partial pressure (temperature). A 
more detail description can be found in the paper by Reuter and Schemer [43]. 
In the limit of low oxygen exposure (partial pressure, chemical potential) on 
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the left hand side of the plot the clean silver surface exhibits the lowest grand 
potential and is hence the thermodynamically most stable structure. When 
the chemical potential increases to around —0.95 eV the thermodynamical 
equilibrium moves to the lowest investigated oxygen coverage (6 = | ML). 
With increasing chemical potential the p(2 x 2) structure at a quarter cover- 
age becomes stable, followed by the two very similar half coverage structures 
(c(2 x 2) and missing row reconstruction). The calculated chemical potential 
of oxygen in the Ag 2 silver oxide structure (i.e. the heat of formation of the 
oxide) of —0.39 eV (experimental value: —0.32 eV [22]) represents the oxygen- 
rich limit for the investigated structures. In other words: for oxygen chemical 
potentials above —0.39 eV it is (from a thermodynamical point of view) more 
favorable to convert the silver surface into Ag20, than increasing the coverage 
at the surface beyond 0.5 ML. 

The result of a maximum oxygen coverage of around 0.5 ML is consistent with 
the experimental estimate of a maximal coverage of around 0.4 ML [3]. Higher 
oxygen exposure leads already to an oxidation of the silver surface. This is 
also in line with the calculated core level binding energies. 



4-2 Core-level energies 

In section 3.8.1 we demonstrated that there is a linear dependence between 
the Ag core level binding energy and the number of oxygen neighbors. The 
comparison with the experimentally determined core- level shift of ~0.4 eV [3] 
implies that most silver surface atoms have only one oxygen neighbor. Only if 
also reconstructions, which reduce the coordination of the surface silver atoms 
are taken into account also two oxygen neighbors become possible (comp. fig. 
5). The results for the oxygen Is binding energies are less simple to interpret. 

Figure 6 summarizes once more the experimental results. The main (low tem- 
perature) species is the hollow adsorbed oxygen species with a coverage less 
than 0.5 ML, which marks also the reference energy of the results compiled 
in table 8. The molecular species is characterized by a higher binding energy; 
Ai^cL =1.4 eV in our calculation compared to 1.7 eV in Ref. [3]. This dis- 
crepancy might be attributed to coverage effects for both the molecular and 
the atomic adsorbates. With increasing temperature Ag vacancies form and 
parts on the surface will undergo an oxygen-induced missing row reconstruc- 
tion; perhaps with additional oxygen in residing in subsurface sites. Increasing 
the temperature even more the mobility within the surface layer also increases 
and the (experimentally observed) phase transition sets in, leading to two new 
peaks in the XPS signals at around 528.3 eV and 530.9 eV, i.e. at —2.0 eV 
and +0.6 eV with respect to the low temperature peak. The peak at —2.0 
eV could indicate low-coordinated oxygen species (bridge, on-top sites). How- 
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ever, the energy difference between hollow and top sites is so large, that an 
occupation of on-top sites is highly unlikely at ambient temperatures. Bridge 
sites represent the saddle points for a diffusion between hollows - the calculated 
core-level shifts, however, are too small to account for the observed low-energy 
peak around 528.3 eV. It is more realistic to assume that the phase transition 
finally leads to the formation of an (oxygen terminated) oxide layer on top 
of the silver surface, characterized by oxygen atoms in bridge and tetrahedral 
interstitial sites of the silver substrate. The calculated values of —1.45 eV and 
+0.57 eV for the perfect silver oxide layer agree reasonably well with exper- 
imental values (—2.0 and +0.6 eV), considering that the newly formed oxide 
is never perfectly flat and will not only consist of (100) facets, but also other 
surface orientations. Furthermore, the lattice constant will be determined by 
a compromise between that of the underlying silver substrate a cxp =4.08A 
(a GGA =4.16A) and that of Ag 2 a cxp =4.72A (a GG A=4.83A). 

4-3 Work-function and frequencies 

The results for work function and stretch frequencies provide further infor- 
mation. However, we are unable to explain the negative work-function change 
observed at low coverage and temperatures in Ref. [1]. The pronounced in- 
crease of the work-function with temperature and coverage is consistent with 
the higher values for A<£> for the lower coordinated sites. 

The adsorbate-substrate stretch frequencies for this particular system also do 
not permit a unique identification of the adsorption geometry. Both, the ex- 
perimental and the calculated values, are biased by high uncertainties. On 
the experimental side the unknown coverage and different preparation meth- 
ods lead to values differing by about 10 %. For the calculation the strong 
coupling between soft surface modes of the substrate with the low vibrational 
frequencies of a relatively heavy adsorbate complicate the analysis. This mode 
coupling could only be taken into account by performing full ab-initio molec- 
ular dynamics simulations over long time intervals. 

4-4 Surface mobility 

The high mobility of the silver surface is an important ingredient for its com- 
plicated behavior. Hence we want to spend a few words on this topic. During 
our study we performed additionally several ab-initio molecular dynamics sim- 
ulations, mainly for simulated annealing optimizations for various structures 
and surface stoichiometries. E.g. by a simulation within a p(2y / 2 x 2) cell with 
a surface stoichiometry of Ag:0=6:3, we were able to reproduce the missing 
row reconstruction. An interesting byproduct of these simulations is the very 
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high mobility of the surface silver atoms, induced by the strong weakening of 
Ag-Ag bonds due to the formation of strong Ag-0 bonds. Although we can- 
not give a quantitative description of this behavior, in all our simulations the 
mobility of the silver atoms was much higher compared to that of the oxygen 
atoms, which were hardly leaving their hollow sites. This observation is further 
supported by STM observations of a high mobility of Ag(100) surfaces under 
atmospheric conditions at room temperature [44]. This mobility favors the for- 
mation of vacancies and other defects well below the melting temperature , 
so that the effective silver coordination of the oxygen atoms is reduced. The 
consequence is an early onset of the metal-oxide transition accompanied by 
a variety of defects and small facettes exhibiting different orientations of the 
metal as well as the oxide surface. All this complicates the comparison with 
well defined surface structures as investigated in a theoretical study. 



5 Conclusion 



We have presented detailed ab-initio density-functional studies of the atomic 
structure, energetics, and other properties of the clean and oxygen covered 
Ag(OOl) surface. Besides adsorption in high symmetry positions also molecu- 
lar adsorption, an oxygen- induced missing-row reconstruction and the Ag 2 
silver oxide surface was investigated. In comparison with data taken from the 
experimental literature we find the following scenario: At low temperature 
oxygen adsorbs molecularly, forming a c(2 x 4) over layer of molecules residing 
in hollow sites. At about 200K these molecules dissociate. The atoms adsorb 
at hollow sites with coverages never exceeding 0.5 ML. Due to the strong Ag- 
O bond, the intermetallic bonds are weakened and partially broken, so that 
vacancy structures such as the missing row structure become possible. At even 
higher temperatures (around 320 K), all atoms at the surface become mobile, 
so that during silver diffusion the coordination number of the oxygen adsorp- 
tion sites is reduced to 3 or even 2. With increasing oxygen exposure and 
temperature the silver surface turns into an Ag 2 surface covered with bridg- 
ing oxygen atoms. By comparison with experimental findings[3], we propose 
that this species is finally the reactive one for e.g. CO or C2H4 oxidation. 
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